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We present a fully nonlinear study of long wavelength cosmological perturbations within the 
framework of the projectable Hofava-Lifshitz gravity, coupled to a single scalar field. Adopting the 
gradient expansion technique, we explicitly integrate the dynamical equations up to any order of 
£N) , the expansion, then restrict the integration constants by imposing the momentum constraint. While 

7—i ■ the gradient expansion relies on the long wavelength approximation, amplitudes of perturbations 

do not have to be small. When the A — » 1 limit is taken, the obtained nonlinear solutions exhibit 
a continuous behavior at any order of the gradient expansion, recovering general relativity in the 
5_l ■ presence of a scalar field and the "dark matter as an integration constant" . This is in sharp contrast 

to the results in the literature based on the "standard" (and naive) perturbative approach where in 
^5 | the same limit, the perturbative expansion of the action breaks down and the scalar graviton mode 

appears to be strongly coupled. We carry out a detailed analysis on the source of these apparent 

■ pathologies and determine that they originate from an improper application of the perturbative 
CN| ' approximation in the momentum constraint. We also show that there is a new branch of solutions, 

valid in the regime where |A — 1| is smaller than the order of perturbations. In the limit A — >• 1, this 
new branch allows the theory to be continuously connected to general relativity, with an effective 
1 component which acts like pressureless fluid. 

& 

I. INTRODUCTION 

<N; 

Recently, Hofava [l| proposed a new theory of quantum gravity in the framework of quantum field theory. One 
of the essential ingredients of the theory is inclusion of higher-dimensional operators, so that they dominate the 
ultraviolet (UV) behavior and render the theory power counting renormalizable. Improvement of the UV behavior by 

■ higher-dimensional operators has been known for some time |2j but in those previous attempts, higher time derivative 
terms led to ghost degrees of freedom. The major modification put forward by Hofava's theory is that the power- 
counting renormalizability is achieved without inclusion of higher time derivative terms. This is realized by invoking 
the anisotropic scaling between time and space, 



ON 

o 



t^b- z t, x^b~ l x, (1) 

so that higher-dimensional operators include spatial derivatives only. This is reminiscent of Lifshitz scalars Q in 
condensed matter physics, hence the theory is often referred to as the Hofava-Lifshitz (HL) gravity. For the 3 + 1 
dimensional theory to be power-counting renormalizable, the dynamical critical exponent z has to be larger than or 
equal to 3 [l[ (see also Q). Because of the anisotropic scaling, the theory cannot be invariant under the spacetime 
diffcomorphism, x^ — > x l>l {x v ), (/i, v = 0, 1, 2, 3). Instead, the fundamental symmetry of the theory is the invariance 
under the so called foliation-preserving diffeomorphism, 

t^t'{t), x^x'(t,x), (2) 

denoted usually by Diff(M, J 7 ). The basic variables of the theory are the lapse function N, the shift vector N l , and the 
3-dimensional spatial metric gij [5j . Since the lapse function N corresponds to a gauge degree of freedom associated 
with the space-independent time reparametrization, it is natural to restrict the lapse function to be independent of 
the spatial coordinates: 

N = N(t) . (3) 
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This condition, imposed in the original formulation of the theory, is called the projectability condition. 

Since its introduction, there has been many cosmological applications of the HL gravity and various remarkable 
features have been found (see jf| 0] for reviews). In particular, the higher-order spatial curvature terms can give rise 
to a bouncing universe Q , may ameliorate the flatness problem Q and lead to caustic avoidance (lpj ; the anisotropic 
scaling provides a solution to the horizon problem and generation of scale-invariant perturbations without inflation 
[ill ] , a new mechanism for generation of primordial magnetic seed field [13] , and also a modification of the spectrum of 
gravitational wave background via a peculiar scaling of radiation energy density [l3j ; with the projectability condition, 
the lack of a local Hamiltonian constraint leads to "dark matter as an integration constant" [14| ; in the parity- violating 
version of the theory, circularly polarized gravitational waves can also be generated in the early universe |15| . 

Despite all of its remarkable features, the theory has been challenged by significant questions. In particular, the 
Diff(M, J 7 ) symmetry allows the existence of an additional spin-0 degree of freedom, often called scalar graviton, and 
its fate is one of important open issues. Actually, the scalar graviton is known to be unstable either in the UV due 
to ghost instability or in the infrared (IR) due to gradient instability [H|-[l8j], depending on the value of a coupling 
constant A. In order to avoid the ghost instability, A must satisfy either A<l/3orA>l. Precisely in these two 
ranges, the scalar graviton exhibits gradient instability at long distances. We then have to tame this IR instability by 
expansion of the universe [19l . |20| or have to hide it by the standard Jeans instability. One can formulate a condition 
under which one of these happens . Essentially, the condition says that A must be sufficiently close to 1 in the IR. 

However, in the limit A — > 1, the scalar graviton appears to be strongly coupled [20l - |22l |. That is, the "standard" 
(and naive) perturbative expansion breaks down in the sense that nonlinear terms dominate linear terms in the A — > 1 
limit. Note that this does not necessarily imply loss of predictability: if the theory is renormalizable, all coefficients 
of infinite number of nonlinear terms can be written in terms of finite parameters in the action, as several well-known 
theories with strong coupling (e.g. .23]) indicate. However, because of the breakdown of the (naive) perturbative 
expansion, we need to employ nonperturbative methods to analyze the fate of the scalar graviton in the limit A — > 1 . 
Such an analysis was performed in Q for spherically symmetric, static, vacuum configurations and it was shown 
that the limit is continuously connected to general relativity (GR). 1 This may be considered as an analogue of the 
Vainshtein effect [24], HH . A similar consideration for cosmology was given in (2(| , where a fully nonlinear analysis of 
super horizon cosmological perturbations was carried out. 

One of limitations of the analysis in [26| is that it is for a purely gravitational system in the absence of ordinary 
matter (but with "dark matter as an integration constant"). Since the naive perturbative expansion is known to 
break down not only in the gravity sector but also in the matter sector [20l [2lj , it is rather important to extend the 
analysis of [26] to the system with ordinary matter. Technically speaking, however, this kind of extention is indeed a 
nontrivial challenge since the system now has multi components (ordinary matter and "dark matter as an integration 
constant") and the gradient expansion technique has not been developed for multi-component systems even in the 
standard cosmology in GR. 

Thus, one of the main objectives of the present paper is to extend the analysis of [26] to the case where HL gravity 
is coupled to a single scalar field, and provide yet another example indicating that general relativity (plus "dark 
matter as an integration constant") is restored in the A — > 1 limit by nonlinear dynamics. Another goal is to point 
out the source of the discrepancy between perturbative and nonperturbative results. As we will see, the solution of 
the momentum constraint in the naive application of the "standard" perturbative expansion is not valid in the regime 
where |A — 1| is smaller than the order of perturbations. 

The paper is organized as follows. In Sec. |Hl we briefly review the basic equations in the HL gravity with the 
projectability condition while in Sec. lIIIl we analyze the inhomogeneous cosmology in HL gravity using the gradient 
expansion method [27| . and present the solutions to the equations of motion. In Sec lIVl we present a discussion on the 
source of the divergences in the naive perturbative expansion and show that the momentum constraint is dominated 
by nonlinear terms in the A — > 1 limit. The results are summarized and discussed in Sec|Vl The paper is supplemented 
by two Appendices, in which we present some of the technical steps of our calculations. 



1 Specifically, the solutions are continuously connected to the A = 1 theory, whose action has the exact same form as the Einstein-Hilbcrt 
term (up to high curvature terms negligible at low energies). However, due to the different symmetries, the resulting theory is not 
exactly GR, but GR with an effective component which acts like dark-matter [141 . This is what we mean by "continuity with GR" 
throughout the present paper. In the case considered in [f|, however, the "dark matter" component is automatically set to zero by the 
assumed staticity. 
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II. BASIC EQUATIONS 

In this section, we review the basic equations of the HL gravity coupled with a scalar field [t| [H[ , following the 
notation in and reformulate them in a way suitable for gradient expansion [26]. In order to make the present 
paper self-contained, some repetition of the material in (2(| is inevitable in SecslIIl and lllll although we shall try our 
best to limit them to a minimum. 

With Diff(M, J 7 ) and the projectability, the building blocks of the theory are g^, Ify, Di and where K^ 
denotes the extrinsic curvature of constant time hypersurfaces, Di is the covariant derivative compatible with the the 
3-dimensional spatial metric , and Rij is the three-dimensional Ricci tensor built out of gij . (This is in contrast to 
GR or any other theory with general covariance whose building blocks are the 4-dimensional metric and its Riemann 
tensor.) For the critical exponent z = 3, their momentum dimensions are, respectively, [Kij] = [k] 3 and [Rij] — [k] 2 . 
Throughout the present paper, we shall impose the projectability condition as well as invariance under the spatial 
parity (x 1 — > —x 1 ) and the time reflection (t —> —t). The number of independent coupling constants in this setup is 11 
for z — 3 [UQi!]. In fact, with the foliation-preserving diffeomorphisms ([2]), the projectability condition (j3)), and the 
additional requirements of parity and time reflection symmetry, the most general gravitational action can be specified 
as 

h = ^ L I Ndt^d 3 x (KijKV -\K 2 -2k + R + L z>1 ) , (4) 
where g is the determinant of , and the extrinsic curvature Kij is defined as 

K * = ~ D i N 3 ~ D 3 N i) , (5) 

K (= g lJ Kij) is the trace of K^, and R is the Ricci scalar constructed from g^. To lower and raise an index, g^ and 
its inverse g n are used. For the sake of simplicity and clarity, in the remainder of this paper, we choose our units 
such that Mpi = 1. 

In contrast to GR, the less restricting symmetry allows both kinetic terms KijK 1 ^ and K 2 to be invariant inde- 
pendently, giving rise to the extra parameter A, which assumes the value 1 in GR, as mentioned above. Furthermore, 
in order to realize the power-counting renormalizability, the higher curvature Lagrangian L z> i should include up to 
sixth spatial derivatives. For the analysis in the present paper, the concrete form of L z> \ is not needed. Adding 
the scalar field action 1$ that is invariant under spatial parity and time reflection, as well as the foliation preserving 
diffeomorphism, the total action is 



i = h + h 



U = J Ndt^d 3 



i(a ± 0) 2 -y(0,A, 5y ) 



(6) 



where we define the derivative along vector normal to the hypersurface 



d ± ^±(d t ~N k d k ), (7) 

and we decompose the scalar field potential as 

V(<f>, D h 9ij ) = V a (0) + V z > x (0, A , 9ij ) . (8) 

Here, V z >i summarizes terms with two or more spatial derivatives and like L z> \ above, its concrete form is not needed 
for the purposes of the present paper. 

Variation of the total action with respect to the 3-dimensional metric g^ leads to the dynamical equation 

£gij + £<f>ij = , (9) 
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where 



9ik9jl 



2 SI g 

N y/9 S 9kl 



- N k D k ) Pt3 + ^( Plk D 3 N k + Pjk DiN k ) 
-Kpij + 2K k p kj + ^g l0 K kl p u - Agij - Gy + £ g , z >i, 



£<f>ij — 



9ik9jl 



9>j 



■,z>l,ij 



(10) 



(ii) 



Here, £ g , z >i,ij and S^. z >\,ij are contributions from L z> i and — V z >i, respectively, pij = Kij — XKgij, and G,*j is 
Einstein tensor of g^. The trace part and traceless part of Eq.© are, respectively, 



(3A-1) [8±K+-K 2 



and 



where 



dxA^ + KA\ + ~^{A k j d k N i - A\d s N k ) 



(d ± <j>) 2 + z = o , 



-zA 

3 3 



o, 



A l ^K^--K5) 



is the traceless part of K l ' ■ and we defined 



(12) 



(13) 



(14) 



Z 9J + 



yi — 



Z = Z\, 
—A5j - G l j + g lk £ z >i, g ,kj 
-V (4>)S i j +g ik £,>i,*,k3- 



(15) 



Here, Z % g • is the variation of the potential part of the gravitational action with respect to the spatial metric; it is 
a generalization of (minus) the Einstein tensor of gtj to include higher curvature terms, as well as the cosmological 
constant. The quantity Z^ • is obtained similarly from the potential part of the scalar field action. 
The variation of the total action with respect to <fi yields the remaining dynamical equation 







1 SI, 



NV9 



dt(y/gd±<j)) 



N V9 



Di(^gN l d A 



Eft, 



where 



1 5 

V9 H 



^/g~dtd 3 x V{cf>,D u g l3 ) = V£(0)+-Ek,,>i 



(16) 



(17) 



and £ , j2 >i is the contribution from V z >\. 

Since the 3-dimensional spatial diffeomorphism is a subgroup of the foliation preserving diffeomorphism, Z l g • and 
Z 1 ^ . satisfy the generalized Bianchi identity and matter conservation, 



0, D j Zl i + Etd i (t> = 0. 
For convenience, we decompose the spatial metric and the extrinsic curvature as 



= a 2 (t) 



1 



K(t,x)S%+A l Jt,x). 



(18) 

(19) 
(20) 
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where we have defined C{t,x) so that det7 = 1, and a(t) (up to an overall normalization) is defined later in Eq. (|37l) . 
The trace part and the traceless part of the definition of the extrinsic curvature lead, respectively, to 

d±C+^ = l (k + Iw*) , (21) 



JVa 3 1 N 

and 



0j.7y = 2 UkA^ + - [ 1]k d t N K + jikdjN" - - jijdkN" ) . (22) 



N y u~ 3 

The momentum constraint is obtained by varying the action with respect to N l : 

D j K j i -\d i K = d ± <f>d i <f>, (23) 
According to the decomposition (|20[) . the momentum constraint is rewritten as 

B 3 A\ + 3 A\ djC - \a\ ( 7 ~ 1 ) lk d i7jk - ^ (3A - 1) diK = drf . (24) 

It can be shown that the evolution equations we have derived are consistent with vanishing A l i , lndet7, 7^ — 7^ 
and ".,/,. 1 A , V-, HJ. 



III. GRADIENT EXPANSION 



In this section, we analyze the dynamics of nonlinear superhorizon perturbations in the spatial gradient expansion 
approach. This approach is valid as long as the characteristic length scale L of the perturbations is much larger than 
the Hubble length By the introduction of small parameter e ~ 1/(H L), we perform a series expansion on all 

relevant quantities and equations. For instance, a spatial derivative acting on a quantity at order e p raises the order 
to e p+1 and thus is counted as 0(e). We then solve the equations order by order in gradient expansion, extending the 
calculations of 26] in a spatially flat Friedmann-Robertson- Walker background, to include a single scalar field as the 
source. 



A. Gauge fixing 

The foliation preserving diffeomorphism invariance, like all other gauge symmetries, reflects a redundancy in the 
descriptions of the theory. By an appropriate choice of gauge conditions, these degrees can be eliminated and physical 
quantities can be extracted. In the present paper we adopt the synchronous gauge, or the Gaussian normal coordinate 
system, by setting the lapse function to unity and the shift vector to zero: 

N = 1 , N l = . (25) 

This choice fixes the time coordinate but in the spatial coordinates, there remains a gauge freedom of time- independent 
spatial diffeomorphism, corresponding to the change of coordinates on the initial constant-time hypersurface. This 
residual gauge degree of freedom will be discussed later in Subsection IIIIF1 

After the gauge fixing, our basic equations (fT2"l) . ([TB")), (TTtJl) . ([2"Tjl and (|2"2"|) are simplified to 

(3X-l)d t K = -l(3\-l)K 2 --A^A^-^) 2 - Z , (26) 

dtA'j = -KA*, + Z l 3 - iz«S* 3 ., (27) 

= d^ + Kdtcp + E^, (28) 

d t { = -^+\K, (29) 
a 3 

d tlrj = 2 7ifc i4*., (30) 
while the momentum constraint (1241) has the form 

<V, + 3 A^-C - \ A^^-^dajk - \ (3A - 1) diK = dttdrf . (31) 

Hereafter, we assume that A 7^ 1/3; this is consistent with the regime of physical interest A > 1, discussed in the 
Introduction section. 
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B. Basic assumptions and order analysis 

We begin by determining the order of all relevant variables. In the limit e — > 0, we expect a universe that looks 
locally like a Friedmann universe, leading to our starting assumption 

dtjij = 0(e) . (32) 

For the scalar field, a similar assumption leads to di<p — 0(e). However, in order to simplify the analysis, we impose 
the stronger condition 

d l <j ) = 0(e 2 ). (33) 
That is, we assume that 4> , which is the leading order term of <fi, is only time dependent: 

^°) = ^°)(t). (34) 
The first assumption (j3"2"| then implies, from Eq. ([3"D|). 

A* j = 0(e), (35) 

leading, using the constraint equation (|3ip . to 

d t K = 0(e 2 ). (36) 
In other words, the zero-th order part of K depends on t only. This fact enables us to define a(t) by 

=tf(°)(=3tf(t)). (37) 

a(t) 

With this definition of a(t), Eq. $29$ leads to 

dtC = 0(e). (38) 
To summarize, the relevant quantities in the analysis are expanded as follows: 

C = C (O) (a ? )+eC (1) ^a ? ) + e 2 C (2) (i,^)+0(e 3 ), (39) 

7« = fi j (Z) + e 7 %\t,x)+e 2 7 g\t,x)+0(e 3 ), (40) 

K = 3H(t) + eK^(t,x) + e 2 K^(t,x)+0(e 3 ), (41) 

A i j = eA (1)l J (t,x)+e 2 A {2) \ J (t,x) + 0(e 3 ), (42) 

= (o) (t) + e0( 1 )(t,f)+e 2 0^(t,f)+O(e 3 ), (43) 

where a quantity with the upper index (n) corresponds to the n-th order term in the gradient expansion. 

C. Equations in each order 

After determining the orders of all physical quantities, we now use this information in the evolution equations 
(126l) - (|30[) to obtain the evolution equations at each order. In the zero-th order of gradient expansion we have 

(3A-1) (d t H + ^H 2 ^ = -i(d t ^) 2 + V Q (^)+A, 
d 2 ^ 0) +3Hd t <f> (0) + V'(<j) {a) ) = 0, (44) 

where a prime denotes the ordinary derivative with respect to the indicated argument. By using the second of the 
above, the first equation can be integrated to give 



9 2 

3H 2 = 



3A- 1 



i(3 t 0(°)) 2 + ^ (O) ) + A 



(45) 
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where C is an integration constant. The last term in the right hand side of this equation is the "dark matter as an 
integration constant" [l4j], a direct consequence of the projectability condition. 
The dynamical equations at order 0(e n ) with n > 1, are written as 



a~ 3 d t 



a 3 ( if (n) + 



3<M n) <9( 



3A- 1 



1 n— 1 



a- 3 8 t (a 3 A {n)l . 



^'(n-p) 

~ p=l 
n-l 



n — l 



2 (3 A - 1) ^ 



^(p)» ^(«-p)j 



TO (0) ) 



3(0 t ^°>) 
3A- 1 



2 1 



3A- 1 ' 

3 j 3 J 



(46) 
(47) 



p=i 



5/ 



3A- 1 



p=i 



3 y 3 A — 1 y 



3A- 1 



(p) 4(«-p) fc 



(48) 

(49) 
(50) 



p=0 



where for later convenience, we introduced new (barred) quantities 

to subtract the terms depending on from (unbarred) zj™^ ! and E^ 1 , defined in Eqs. ljPo)) and (ITT)) . Here, Z^ 1 -, 
Z(") and £?J, n) are the n-th order terms of Z*-, Z, E^, respectively. With this definition, Z^™ ^, Z^™) and do 



nor on 



not depend on C (n) , l\f , K^ n \ A (n) ' p 
. (|3"TT) we obtain tl 

n .. j - 

9^\ + 3 E^ (p) W M - ^E A<P) i^ I ) (,)Iti) i 

p— 1 p— 1 q—0 

4 ( 3 A - 1)3 ,( K « + M^)-g at ^) a 

^ ^ p=i 



Similarly, from Eq. (1311) we obtain the order C(e" +1 ) (n > 1) momentum constraint as 

n n—p 

2 



4(» 



P-<?) 



0, 



(52) 



where (7 1 )(™)^ is the n-th order term of the inverse of 7^, i.e. the inverse (7 1 ) y is expanded as 

( 7 -i)y = /*i + e ( 7 -i)(i)tf + e 2( 7 -i)(2)tf + ... ( (53) 

where = ( 7 - 1 )(°)^ is the inverse of /y. It is straightforward to show that (7 ^("W (n > 1) satisfies the following 
differential equation: 



9t ( 7 -i)(«)tf = -2^A^ i fe ( 7 - 1 ) 
p=i 



(n-p) fej 



(54) 



In addition to the dynamical equations and momentum constraint, there are also some useful identities. First, we 
expand the generalized Bianchi identity (|18p to obtain 



p=i 



/ I \ 1 " n_ P 

d 3 z [n) \ + 3 E (z^\ - i z^ a,c^) - i E E ^fr- 1 )^*^ 



-9) 



D =l 9 =0 



n — 1 

E [e [ ;~ p) + k,"(0 (o) ) ^ ( "~ p) ] 9^ = , 



(55) 



p=i 
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for n > 1. Next, expanding the conditions A l i = 0, ftlndet7 = 0, 7^ — 7^ = 0, likA k j — r fjkA k i — and 
A l j — 7jfe^4 fe ;(7~ 1 ) /i = leads to the following identities: 

n 

A {n) \ = , ]£(7- 1 ) WJfc *#- p) = , 7? - 7? = , 

n— 1 n— 1 n— p— 1 

E^'^r^'^C ^"VE E T^™^^ 1 )^^. (56) 

p=0 p=0 q=0 



D. 



(e) solution 



For 0(e), Eqs.(gBJ(-J5DJ| reduce to 

9i 



3 (>> , 3£V^ 
1 3A-1 / 



0. 
0. 



3 a t ( a 3 a t <?!) (1) 



TO (0) ) 



3(9 t 0(°))- 
3A- 1 



0« 



* 3 \ + 3A-1 J 3A-1 ' 



8?$ = 2/ lfe A 



(i)k 



(57) 
(58) 
(59) 

(60) 
(61) 



where from equations (IT51) . (IT7|) and (|51l) . we have at first order, • — = = 0. Integrating the above 



equations, we obtain 



K 0) 



3 0«c^°) CW(x) 



3A- 1 



A 



a(t) 3 

i 

j v 



0(t)3 ' 



a{t'fW{t') 



A ln a(t') 3 ^(i') +0inW 



C (D 



7 (1) 



dtf 



1 



2/ ifc (f)c (1)A ;.(f) 



a 3 (i') 3A-1 



AW 
AW, 



3/^/\ + 7inij C 3 -) 



(62) 
(63) 

(64) 
(65) 
(66) 



where the integration "constants" , C^ 1 ^, <p^ , </>[„', an( i 7ini/ depend only on the spatial coordinates af 1 



and satisfy 



C" 



(i), 



(i)k 



(67) 



The functions /j(t) (i = 1, 2) are two independent solutions of the homogeneous equation 



a- 3 d t (a 3 d t fi) + 
and 



TO (0) ) 



3(a t </»(°)) 2 



3A- 1 



/i = Q; AWn) = l, AWn)=0; AWn)=0, A(t in ) -1, (68) 



W(t) = h{t)d t f 2 {t)-f 2 {t)d t h(t). 



(69) 
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The two first order integration "constants", and 7^^, can be absorbed into their zero-th order counterparts, £^ 
and 7jny. Thus, without loss of generality, we can set 



c (1) = 



7 (1) .=0. 

'in 11 



(70) 



Finally, the momentum constraint equation (|52p with n — 1 leads to the following relation among the remaining 
integration constants, C^ 8 ., and /y, 



+ 3C (1) VK(°> - ^ c (1) \f lk dj Jk - ± (3 A - 1) 9,C« = 0. 



(71) 



Note that <^;j^(x) and $1'' (x) do not appear in this equation. The physical meaning of (j)^' (x) and (/}\^'(x) are obvious: 



AW 



(72) 



E. 0{e n ) solution (n > 1) 



Equipped with the zero-th and first order solution, we can now determine the general solutions at arbitrary order 
in gradient expansion. For any n > 1, the solution to Eqs. (|46)) -([50 |) is 



K (n) + 



3 d t (j)^ 
3A- 1 



3(t) 



n-l 



.4 



dt'a 3 (t') 



£ [A (p) V^^^A ( "- p) fK^',f) + ^^)(^^f)a t ^("-rf(^^x)] , (73) 



n — 1 

y ^(t',f)A ( "" p) \.(t', a -) + z ( " ) ;.(t',f) - -z^{t , ,x)s i 



= 1 



W(f) 



W(t') 



*(«) 



» 



dt' 

3 



3A- 1 



3A- 1 



•/tin p=0 



,(74) 

(75) 
(76) 
(77) 



where 



r (n) (t, 3) - ( * <»> (*, 3) + ) ft^> (t) + £ /fW (t, f ) ft^Cn-P) (t , f ) + (t> f ) , (78) 



and by redefining C^, C*' fc', and fij, we have set, respectively, 



K (n) 



f=ti, 



5 t <?!> (Tl) = C (n) 



(n) 
lij 



(79) 



We remind that the first order constants have already been fixed in Eq. (|70p by redefinition of Q ' and /«, respectively. 



The initial condition for (n > 1) implies that 7y| t=t . - / , , 
(n > 1). Therefore, for n > 1, the solution to Eq. (l5H) is 



/«, (7 _1 ) ij 'L_,. = / ij and ( 7 - 1 )(«)> 



-l\(n) ij 



2/ dt'Y,A (p)i .h- 1 : 

• 7 *.n P =l 



(ra-p) fej 







(80) 



As shown in Appendix IA1 the solution (j?3")) -(fTT |) automatically satisfies the (n + l)-th order momentum constraint 
equation (|52")l , provided that the redefined integration constants (C^ 1 ^, C' 1 fij) satisfy ([7T|) up to C(e™ +1 ). 
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F. Number of physical degrees of freedom 



The solution we obtained in the previous subsection involves a number of functions depending only on spatial 
coordinates, C (0) (^)> C«(£), C (1) V(f), ^(x) and ^\x) which emerged as integration "constants" . 

However, not all of the components are independent nor physical. Firstly, they are subject to the constraint (|71[) . 
Secondly, as stated just after Eq. (j2"5)) . our gauge condition (|2~5|) leaves time- independent spatial diffeomorphism as a 
residual gauge freedom. Therefore, the number of physical degrees of freedom included in each integration "constant" 
is 



C<°>(f) 
C (1) V(f) 



1 scalar growing mode = 1 component , 

2 tensor growing modes = 5 components — 3 gauge , 

1 scalar decaying mode = 1 component , 

2 tensor decaying modes = 5 components — 3 constraints , 
2 scalar modes . 



(81) 



This is consistent with the fact that the HL gravity includes not only a tensor graviton (2 propagating degrees of 
freedom) but also a scalar graviton (1 propagating degree of freedom) and that our system includes a scalar field (1 
propagating degree of freedom) as well. 



IV. PERTURB ATIVE VS NONPERTURBATIVE APPROACHES 



In the previous section, we have derived solutions for nonlinear perturbations in any order of gradient expansion. 
While gradient expansion relies on the long wavelength approximation, amplitudes of perturbations do not have 
to be small. Thus, our analysis in the previous section is totally nonperturbative with respect to amplitudes of 
perturbations. The dynamical equations and their solutions do not suffer from any divergences in the A — > 1 limit, 
and GR coupled with a scalar field and dark matter is safely recovered in this limit. 

This is in sharp contrast with results known in the literature based on the "standard" (and naive) perturbative 
approach, in which pathologies such as divergences and strong coupling are found in the A — > 1 limit. In this section, 
we shall see how this problem arises in the "standard" perturbative approach and why it becomes under control in 
our nonperturbative approach. In the "standard" perturbative approach, all relevant equations are expanded with 
respect to amplitudes of perturbations, irrespective of sizes of coefficients in the expansion. We shall see that, in 
the momentum constraint, coefficients of terms linear in perturbations actually vanish in the A — > 1 limit and thus, 
for sufficiently small but nonvanishing |A — 1|, the linear terms become less important than nonlinear terms. Hence, 
neglecting nonlinear terms, blindly solving the linearized momentum constraint and then taking the A — > 1 limit 
would be totally nonsense and lead to inconsistencies. This is precisely the situation in the "standard" perturbative 
approach. Clearly, this is a breakdown of the treatment based on the "standard" perturbative expansion but not 
of the HL theory itself. Indeed, as already stated above, our nonperturbative analysis in the previous section does 
not show any pathologies in the A — > 1 limit. In the rest of this section, we shall investigate these issues explicitly. 
For simplicity we shall consider the cases without the scalar field (but with the built-in "dark matter as integration 
constant"). 



A. Breakdown of standard perturbative expansion in the limit 

In this subsection let us briefly review the standard perturbative approach and see that, contrary to the nonper- 
turbative approach based on the gradient expansion in the previous section, it breaks down in the A — > 1 limit. 
Let us adopt the following metric ansatz in the transverse gauge, 

2V=1, N^&B + n,, 9ij = a 2 e 2iT (e h ).. , (82) 

where n, is transverse and hij is transverse and traceless: d % rii = 0, d % hij = and K l i = 0. Throughout this 
subsection, indices are raised and lowered by and Sij. We introduce a small parameter e, consider £x, B, rii and 
as quantities of 0(e), and perform perturbative expansion with respect to e. 
In the regime of validity of the standard perturbative expansion, in order to calculate the action up to cubic order, 
it suffices to solve the momentum constraint up to the first order, which can be written in the form, 

di [a 2 (3A - l)d t C T - (A - 1)AB] + ^An* = 0, (83) 
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leading to 



where A = d l dj 



3A — 1 

a- 2 AB = -j—^dtCr, «j = 0, (84) 



(2) 

It is straightforward to calculate the kinetic action up to the third order. The quadratic part IL' and the cubic 



P art 4fn are 12£ 



1 



k 



in 



7, (3) = / dtd 3 xa 3 



3Ct ( aT 2 d t (,T AB + \d t h ij d t hi- \ + la~ 4 ( T d l {d t BAB + Z&Bdid-jB) 
8 / 2 



-^{a- 2 d k h ll d k B - 3d t h ij ( T ) a - 2 d l d j B - -^aT 2 d t h 11 d k h lit d k B 



(85) 



When B is eliminated by using (|54"|) , one can easily see that the quadratic part written in terms of (t — y 2 ^L^ (t 

(3) ~ 

is regular. On the other hand, the cubic part I kin written in terms of C,t is divergent in the limit A — > 1. Thus, the 
perturbative expansion breaks down in this limit. More precisely, the regime of validity of the standard perturbativc 
expansion is 

|CT|«min(|A-l|,l), (86) 

and disappears in the A — > 1 limit. 

Evidently, the breakdown of the standard perturbative expansion in the A — > 1 limit originates from the denominator 
A — 1 in the solution (|84[) to the linearized momentum constraint. 

B. Transformation from transverse to synchronous gauge 

In the standard perturbative approach summarized in the previous subsection, we have adopted the transverse 
gauge (I82[) . Instead, in the nonperturbative approach based on the gradient expansion presented in Sec lIIIl we have 
adopted the synchronous gauge (|25[) . In this subsection, we shall investigate the spatial coordinate transformation 
between the two gauges. (Note that in both gauges the space- independent time reparametrization is already fixed by 
the condition N = 1.) The transformation is nonlinear but we treat it perturbatively. As we shall see below, this 
provides an alternative way to sec the breakdown of the standard perturbative expansion. 

As described in Appendix [Bj we start with the transverse gauge, carry out the spatial gauge transformation to the 
synchronous gauge, and use the momentum constraint (in the transverse gauge) to eliminate the nondynamical degree 
of freedom. In this way, we can express the perturbation in the synchronous gauge in terms of that in the transverse 
gauge. Up to the second order, the result is 

c = -fx^( CT -f^ (9iCT)(9iA " 1CT) + |^iy f dt ' a-^ct) 



(3A-1) r* 



2(A-1) 

dt'A- 1 



dt'A- 



2 (c^ACtX^A-^Kt) + ( d l d^ T + - Aft« ) (c^A-^Cr) + (A^X^'Cr) 



+ 0(e 3 )|, (87) 

where £ is the perturbation in the synchronous gauge defined in (|19j). Ct and hij are the metric perturbations in the 
transverse gauge defined in (|82l) . One can easily see that the terms quadratic in £r are suppressed with respect to 
the linear term under the condition (|86|) . 

Conversely, for a fixed amplitude of £<r, the expansion with respect to e in (|87[) breaks down in the A — ¥ 1 limit. 
This is very similar to the way how the standard perturbative expansion of the action (j85[) breaks down in the A — > 1 
limit. It is apparent from the intermediate steps of the calculation (shown explicitly in Appendix [Bj that the terms 
with negative powers of (A — 1) arc introduced by the solution of the momentum constraint. 
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C. Linear vs nonlinear terms in the momentum constraint 



Having understood that the origin of the breakdown of the standard perturbative expansion is the treatment of the 
momentum constraint, we now discuss the regime of validity of the standard perturbative expansion in the momentum 
constraint. Importantly, we shall see that a new branch of solution emerges at the edge of the regime of validity of 
the standard perturbative expansion. 

For this purpose, we adopt the transverse gauge (|82[) and expand the momentum constraint with respect to C,t and 
hij , considering them as small quantities but keeping B and rii as nonlinear quantities 2 : 

Cr = 0(g), h lJ= 0(q), B = 0(q°), n i = O{q ), (88) 

where we have introduced a small parameter q to count the order of perturbations Qt and hij . In the absence of the 
scalar held, the momentum constraint is [261 ] 

= Ui = djA*. + 3 A\ d 3 ( T - \A\ ( 7 ~7 fc d iljk - \ (3A - 1) diK, (89) 

where 7y = (e h )ij, (7 _1 ) u = (e - ' 1 ) 1 -' , while the trace of the extrinsic curvature (|2ip becomes 

^ = 3(^ + ^)->[(^'iV,] 

= 3H- a- 2 AB + 3cKt + a' 2 [-(d k ( T ){d k B) + 2( T AB + h kl d k diB] 

+0(q) x n k + 0(q 2 ), (90) 



and the traceless part (f22|) is 

A\ = lh-y k d^ kt - ± {(-y-y^udk [(g-^N^ +d t [{g-yNk] - \s{d k [{g- x ) hl Ni] 
= \d t h 3 % - ± fyd j m + dm*) + (tfdiB - ±6{AB 

(d j ( T )(diB) + (diCrX^B) + 2C T (d j d i B) - Isj [(d k ( T )(d k B) + ( T AB] 



(dih jk + 0V - d k h^ (d k B) + h jk (did k B) - yjh^id^B)^ 

+0(q)xn k + 0(q 2 ). (91) 



i ri 



Here, it is understood that (g 1 ) u is the inverse of gij, that derivatives do not act beyond parentheses and that indices 
are raised and lowered by S 1 ^ and <5y . A straightforward calculation results in the following expansion of Hi, 



H, = -{Z\-l)did t (T + 0{q 2 )- ^ [A + 0(q)]m 

2a z 

+ h { (A _ 1} [ s * A + ° {q) } + (I AhJ * + 91 Ut + 6i A(t ) + ° {q2) ] djB ■ (92) 

Notice that in the above, no assumption is made for B and rii, which are still considered to be nonlinear quantities. 
It is now clear that the leading term in the coefficient of B relies not only on the order of perturbations, but also on 
the value of A — 1. 

In the regime q <C min (1, A — 1), the momentum constraint becomes 

O \ _ 1 

a- 2 AB=- -d t ( T + 0(q 2 ), ni = 0(q 2 ), for q < min (1, A - 1) , (93) 

A — 1 



2 Although the discussion in this subsection employs the transverse gauge, the general argument holds in any gauge. In particular, in 
the synchronous gauge, where B = m = 0, the issue arising in B gets transferred to the longitudinal part of hij . However, for the sake 
of clarity, we chose to keep find solutions for the nondynamical fields, gauging away the longitudinal mode instead. 
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and agrees with the result of the standard perturbative expansion (|84)l . Naively using this expression in the action, 
then taking the A — > 1 limit would lead to breakdown of the standard perturbative expansion as already seen in (|85[) 
and (j87| . 

On the other hand, if A is sufficiently close to 1 and the condition 

A - 1 « 5 « 1 (94) 

is met, then the coefficient of B in Eq. (j9"2")) is dominated by the 0(q) terms instead of the 0(X — 1) term. Note that 
this is a nonlinear regime but is still consistent with the assumed smallness of the metric perturbations £r and hij. 
In this regime, the constraint can be written as 

Uj = - + 0(q 2 ) - ^2 (A + O(q)) nj + ^ [Mj. 1 + 0(X - 1) + 0(q 2 )] diB , (95) 



where we have defined 

1 
2 



M/ = - Ah'j + d%C T + 5) ACt = O(q) . (96) 



The transverse part of the Eq. (l^5|) can be found by evaluating d[ k Hj], 

d [k A nj] = (Ah* y )(d k] diB) + (d [k Ah* rfidiB) 

+2(d%( T )(d k] diB) + 2(d [k A{ T )(d j] B) + 0(q 2 ). (97) 

On the other hand, the longitudinal part can be computed from &*Hj as 

-2Ad t ( T + ^MB = 0, (98) 

where we define the operator M as, 

M = M l3 didj + 2 (d l A( T ) $ . (99) 
If this operator is invertible, then we obtain 

B = 2a 2 M- 1 AdtCr + O (— - ) + 0{q) . (100) 



We note that either the (A — l)/q or q term can provide the largest correction to B, depending on the value of A — 1. 

In summary, we have seen that there are two branches of solutions to the momentum constraint, depending on 
the value of A — 1. One is (f9"3"|) in the linear regime, and the other is (|100j) in the nonlinear regime. The standard 
perturbative expansion in the previous subsections corresponds to the solution (|93p . On the other hand, what is 
relevant for the nonperturbative recovery of GR (plus "dark matter") in the A — > 1 limit is the solution f|100[) . The 
two regimes are mutually exclusive. 



D. Yet another consideration 



In the previous subsection we have seen that there are two mutually exclusive branches of solution to the momentum 
constraint. This explains the reason why the standard perturbative approach breaks down in the A — > 1 limit and 
why the theory itself can be still regular and continuous in the limit. 

For A away from 1, the standard perturbative expansion is valid in the transverse gauge and we have the expansion 
of the kinetic action as given in subsection IIV Al On the other hand, for A sufficiently close to 1, i.e. in the regime 
(|M)l . it is the nonlinear solution (|100l) that should be substituted to the kinetic action. 

Unlike the kinetic action, the potential part of the action does not depend on A, when written in terms of £r and 
h^. This is because B does not appear in the potential part of the action. Therefore, if we could somehow define 
a field Cc in such a way that the series of terms in the kinetic action for £t sums up to form a standard canonical 
kinetic term for £ c , then the whole action written in terms of Cc should remain finite in the A — > 1 limit. Since each 
term in the kinetic action (after eliminating B) includes exactly two time derivatives, such a field redefinition should 
be possible in principle. In practice, however, the field redefinition is not easy to perform since it would be nonlinear 
and highly nonlocal in space. Nonetheless, this consideration already suggests that the A — > 1 limit should be regular 
and continuous nonperturbatively. 
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V. SUMMARY AND DISCUSSION 

In this paper, we have performed a fully nonlinear analysis of superhorizon perturbations in the HL gravity coupled 
to a scalar field, by using the gradient expansion technique 27]. After applying the long wavelength expansion to the 
set of field equations, we integrated these explicitly to the second order. We then showed that the solutions can be 
extended to any order in gradient expansion, while satisfying the momentum constraint at each order. These solutions 
are continuous in the GR limit A — > 1 for any order in the expansion, both in the gravity sector, which consists of the 
"dark matter as an integration constant" , and in the matter sector, which contains a scalar field in the present work. 
The form of the equations suggests that our qualitative result should remain the same when additional matter fields 
are introduced. 

This is in sharp contrast with the results obtained in the framework of the "standard" (and naive) perturbation 
theory, in which pathologies such as divergences and strong coupling are found in the A — » 1 limit 0, HH . We 
determined that the results of the standard perturbative expansion are valid only in the region where |A — 1| is larger 
than the order of perturbations. In other words, the range of validity (|86p of these solutions has zero measure in the 
limit A — > 1. We found that the divergences are originating from the momentum constraint, where the coefficients 
of the terms linear in perturbations vanish in the A — > 1 limit. Thus, for sufficiently small but nonvanishing |A — 1|, 
the linear terms become less important in comparison to the nonlinear ones. Neglecting nonlinear terms and naively 
solving the linearized momentum constraint, then taking the A — > 1 limit turns out to be the main source of the 
said pathologies. Once their origins were understood, we carried out a detailed analysis of the nonlinear momentum 
constraint in the perturbative approach. In addition to the known result which is valid when |£| <C min(|A — 1|, 1), we 
found a second branch of solution valid in the regime where |A — 1| <C |£| <C 1). The presence of the latter solution 
justifies the recovery of GR obtained in our nonperturbative approach, in the A — > 1 limit. 

Our results, together with the similar examples studied in [f| [2(| , discernibly support that the apparent strong 
coupling found previously in the HL gravity may only indicate the breakdown of the treatment based on the naive 
perturbative expansion but not of the theory itself. General relativity should be recovered by nonlinear effects, 
analogously to the Vainshtein mechanism [24| that were first encountered in the massive gravity theories. 

We note that the present analysis was limited to the discussion of the classical (and superhorizon) evolution of 
perturbations; their quantum mechanical origin were not considered. In contrast, in Ref.|29[, the HL gravity in the 
A — > oo limit was found to be weakly coupled under a certain condition, and the spectrum of perturbations that were 
generated from quantum fluctuations was calculated in this limit. However, the presence of regular behavior both in 
A — > oo and in A — > 1 limits is not sufficient to draw conclusions on the transition between the two regimes. This is 
because of our lack of an understanding on the details of the renormalization group (RG) flow. Specifically, to be able 
to match these two results, one needs to define a conserved quantity (like the comoving curvature perturbations in 
relativistic cosmology). However, since the matching involves a wide range of varying A, one needs to know how the 
flow of A is realized and how such a flow affects the evolution of cosmological perturbations. With these considerations, 
we refrain from exploring the cosmological implications of our results for now. 

On the other hand, a quantum mechanical extension of our analysis may have a chance to address such issues. 
One of the major concerns with a proper renormalization analysis in HL gravity is the strong coupling problem in 
the A — > 1 limit, or more specifically, the breakdown of the perturbative expansion. However, we have shown in 
Sec. IIV Cl that the full nonlinear analysis in the limit A — > 1 is still consistent with small perturbations, except for the 
nondynamical mode B which becomes nonlinear. The solution to the momentum constraint in the two regimes, (|93[) 
and (jlOOp . gives the nondynamical mode as 

r §j£l ^A^dtC = 0(C) , for |C| < min(A - 1, 1) 
B~\ . (101) 

I 2 a 2 M" 1 Ad t C = 0(1) , for A - 1 < |f| < 1 

Note that both cases are compatible with small (. Thus, substituting this nonlinear solution for B in the action, then 
applying the perturbative expansion for £ may provide a healthy perturbative action. (However, the reduced action 
is nonlocal in space while it is local in time). 

We also note that the breakdown of the naive perturbative expansion does not necessarily result in loss of renor- 
malizability. We know that in the regime |£| <C min(A— 1, 1), the leading UV contributions in the action are invariant 
under the scaling ([1]) with z = 3, provided that the scalar graviton is assigned a vanishing scaling dimension, i.e. 
(, — > £. This fact is nothing but the power-counting renormalizability of the theory, and is seen after replacing B 
in the action with the linear solution (f93|) . or the first line of (|101|) . Note also that coefficients of all possible terms 
in the perturbative expansion are expressed in terms of 11 coupling constants in the action (j4|). On the other hand, 
for the regime A — 1 |^| 1, we need to replace B in the action with the nonlinear solution (jlOOp . or the second 
line of (jlOip . What is important here is that the scaling dimension of B from the nonlinear solution and that from 
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the linear solution are exactly the same: B — >• bB under the scaling ([lj with z — 3 in both cases. Therefore, after 
substituting the nonlinear solution to B in the action, we still conclude that the leading UV contributions in the 
action are invariant under the scaling ([TJ with z = 3, provided that the scalar graviton is assigned a vanishing scaling 
dimension. In other words, the conditions for power-counting renormalizability of the theory continue to hold in the 
nonlinear regime. 

In the present paper, we have considered the projectable version of the HL theory and showed that the general 
relativity (plus dark matter) is safely recovered in the A — > 1 limit. If we relax the projectability condition and thus 
allow the lapse function to depend on spatial coordinates then we should include as the building blocks of the theory 
not only gij, Kij, Di and Rij but also at = c\ \n(N) (with [oj] = [fc] 1 for z = 3) [13] • This gives rise to a proliferation 
of independent coupling constants. For example, for the minimal value of the dynamical critical exponent z = 3, the 
number of independent terms in the gravitational action of the non-projectable extension turns out to be more than 
70 [22f | . In some regime of parameters, the non-projectable extension is claimed to be free from the breakdown of the 
standard perturbative expansion method in the A — > 1 limit, while in other regime the expansion breaks down. It is 
certainly worthwhile performing a nonperturbative analysis of the non-projectable theory in the regime of parameters 
where the standard perturbative expansion breaks down and then identifying the observationally viable regime of 
parameters. 

There is yet another extension of the HL theory, with an additional local (7(1) symmetry, U(l) k Diff(Af, F) 30]. 
It has been shown that the standard perturbative expansion does not break down in the gravitational sector, but does 
break down in the matter sector, at least apparently 31 1 (see also [32T436T ] for more on this extension). It is intriguing 
to see if a nonperturbative analysis similar to those presented in the present paper can resolve this problem. 

Finally, the biggest obstacle in front of the HL gravity, and in general, any Lorentz symmetry breaking theory, is 
the restoration of the Lorentz invariance in the matter sector at low energies [37], HH . Even if the Lorentz violation 
is restricted only to the gravity sector, the radiative corrections from graviton loops will generate Lorentz violation 
in the matter sector. Such terms can be under control provided that the Lorentz breaking scale is much lower than 
the Planck scale [3{| . Another approach is to introduce a mechanism, or symmetry to suppress the Lorentz violating 
operators at low energies, such as supersymmetry 40]. Such an approach was adopted in [4l| where a SUSY theory 
with anisotropic scalin g w as constructed. On the other hand, this seems to be a highly nontrivial task for the case of 
interacting models [421 143| . 
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Appendix A: Order C(e n+1 ) momentum constraint for n > 2 

In this Appendix, we prove by induction that the order C(e") solution ([73")) -([77 ]) satisfies the order 0(e n+1 ) mo- 
mentum constraint equation (|52[) for n > 2. 

The proof extends the method presented in Rcf. 26] to include a scalar field source: we rewrite the left hand side of 
the (n + l)-th order constraint (|52|) as a linear combination of lower order constraints by using the explicit solution 
(IT3"|) - ([77|) . To achieve this, we make use of the generalized Bianchi identity (1551 as well as the identities in Eg. (1561 . 
We also use the following identity for functions f(t) and g(t) satisfying a 3 (i; n ) /(£;„) g(ti n ) = 0, 

f(t)9(t) = J^f t rf *'« 3 (*0 W)- Z dAa^t')f{t'))-g(t') + f{t')-d t ,g{t')\ . (Al) 

By applying the identity (XT) to (/(*),»(<)) = (A (p) \, d^^), (f(t),g(t)) = (A ip) \, ( 7 -l)<*) lk d a %- p ~ q) ) and 
(f(t),g(t)) = (dt(f>( p ' , di(j)( n ~ pS> ), the left hand side of the (n + l)-th order momentum constraint equation ([52]) is 



rewritten as 



n 1 n n—p 



p-q) 



n-1 



p=i 



1 



• 1 „3 



V(*') ^ 3 E [a- 3 ^ (a 3 4 Wi i) d^ n ~^ + A {p)i ^ (d t >C (n - p) 
I p=i 



*M «.„,("-P-9) 



p=l g =0 

n-1 



p=i 



Using Eas.([75 ]) -([7i )) . ([2g ]) -(|50 ]l and ([55]). this is further rewritten as 



,(n+l) 



n / p-1 \ n-1 , 

5] ^^^M 9j ( M +E^' A 



p-i 
E 

.p=2 \ q=l 
n n—p 



p=l 



^EE 

p=l <j=0 



E# (r) ^ (p_r)i i I (7 _1 ) {,)Ifc a i 7j n " p ~' ) 



'j A- 



-2E^ (r)i m (7- 1 ) (9 - r)wfe Wr p ~ 9) 



/ n— p— g— 1 \ 
+A (P)i i(7 -l)( g )« 9 . 2 £ 7 (r) A (n-p-«"0- 



k 7 



i(3A-i)fl, f-E^ (p) ^ p) ) f-E 

V p=i / V p=i 

n — 1 p 1 ™ ™ — p 

E E^ (9)9 ^ (p ~ 9) ^ (Tl ~ p) + - E^ (p) E^ 1 )^^^. 



(n-p-q) 
jk 



p—1 q—1 p—1 q—0 ) 

where we have used the generalized Bianchi identity (|55[) . By using the identities (|56[) we finally obtain 

n-1 

/ J+>n*(+" 

a 3 ' 



c 



(n+1) 



i- / d iV(oE A ^ p) <t +1) - 

w J t ia p = l 
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(A2) 



(A3) 



(A4) 



Since the 0(e 2 ) constraint in Eq. (17T]) is already satisfied, i.e. Q = 0, the above relation implies that C^ n+1 ' = for 
?i > 2. 



Appendix B: Expansion of the nonlinear perturbations 

Here, we present the detail of the calculations to obtain the expression of ( in terms of C,t and , given in Eq. (1571) . 
While the former £ is defined in the synchronous gauge iV,- = 0, the latter £t and hij are defined in in the transverse 
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gauge 5 dkhtj = 0. In both gauges, the freedom in the time coordinate is fixed by the choice N = 1. For the 
perturbative expansion of the spatial metric, we use 

g l3 = a 2 e 2C (e h ) = a 2 6 tJ + I a 2 (2^y + &y) I + 



(AC 2 S i3 +ACh lJ +h a h l J ) 



0(e 3 



(Bl) 



where e denotes the order of perturbations and the indices of hi 3 are raised and lowered with Kronecker delta. 
Throughout this Appendix, when the expansion of a quantity is shown, the terms outside parentheses, in parentheses 
and in square brackets are of order e°, e 1 and e 2 , respectively. 



1. Expansion of the momentum constraint 



We first concentrate on the linear perturbation in the transverse gauge. To remove the nondynamical degrees in 
the shift vector, we solve the constraint equation order by order. We expand the shift vector while separating the 
contributions from each order, as 



T(l) 



N, 



T(2) 



With these decompositions, we expand the momentum constraint in vacuum 

H j = D i K i j -Xd j K = 0, 
as a series in perturbations. At first order, we get 



HP = - (3 A - 1) d.dtCr + ^ fl^ABW - -^AA ? T(1) 
3 a 2 2 a 2 3 



A- 1 



1 



which can be solved by 



ATf (1) = . 



Using the second of these results, the next order constraint yields, 



nf = (a - 1) dj 



- AB ( V + (d k ( T ) (d k B^) - 2 Ct Afif 1 ' - h kl d k diB^ 



[(dth kl )(d k h l3 ) - h kl d k d t h l3 - {d t h kl ){d 3 h kl )} + ^ (d k ( T ) (d t h kj ) 

+ (d k d 3 ( T ) (d k B^) + (ACt) (d 3 B^) + I (Ah k3 )(d k B^) - \ ANf {2) = 



(B2) 

(B3) 
(B4) 

(B5) 



(B6) 



For the following, only the longitudinal part of this relation is needed; by taking its divergence, then using the first 
equation of (|B5p . we end up with 

3A- 1 



AB™ = 2a 2 (^— Y \[2CTd t CT + h^d l d 3 A~ 1 d t CT~(dXT)(d l A- 1 d t C T ) 



2 a 2 (3 A - 1) 



(A -if 

„2 

■A" 1 



A" 1 



2 (^ACtX^A-^Ct) + ( d l d 3 Q T + ^ ) (ftfljA-^Cr) + (ACt) (9*Ct) 



A - 1 



l - (did t h jk ) {d l h 3k ) + l - (d t h l3 ) (Ah 13 ) - 3 (fi^-Cr) (^) 



(B7) 



2. Coordinate transformations 



Next, we determine the transformation between the transverse and synchronous gauges. We parametrize the 
coordinate transformation as 



0(e 3 ), 



(B8) 
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where over-tilde denotes quantities in the synchronous gauge. The parameters M are decomposed as 



(B9) 



with di^ nS> = 0, while the indices of are raised and lowered by S 1 ^ and <5jj. For any tensor field expanded as 



(«r) + 









+ 0(e 3 ), 

the transformation at linear and quadratic order proceeds through [44| 

^ = c5 2 T + 2£ ?(1 ,(5T + £2 (i)r (o) + ^ (2 , r (o). 
For the metric tensor, the transformations become 

Sg^, = 5 9fiV + g$d^ ° + ° , 



(BIO) 



(Bll) 



(B12) 



We now determine the transformation £^ needed to go from the transverse gauge to the synchronous gauge. For this, 
we evaluate the Oi components of (|B12|) and set Sgoi — 0. At first order, we obtain, 



dtBW + N? (1) + a 2 + ) = 



(i) 



where the transverse and longitudinal parts can be easily separated to give, 



dt 



Using the solutions (|B5[) to the linear momentum constraint, the transformation parameters become 



(B13) 



(B14) 



(B15) 



Similarly, the Oi component of the second order transformation (|B12|) gives, 

d l B { ^+Nf {2) W (d 2 cW 2) +dtC- 2) )+{d ] C (1) ) {didjB^-idid^W) (d j BW)-4C T d i BW-2hi j d i BW = 0, (B16) 

where we used the second equation of (|B15|) . Using also the first equation of (|B15[) as well as the first order constraint 
(|B5|) . the longitudinal part of the second order transformation can be obtained as 



AB' 2 » /3A-1 



3A - 1 



3A- 1 
A- 1 



dt' 



h? j a i d j A- 1 d t ,c T 



x A- 1 



(B17) 



Inserting the expression of B^ from (|B7I) into the above expression, we finally obtain 



A £(2) = 



3A- 1 



4(3A-1) 



— {dXr) (diA-'Cx) - ( \_ 1)2 ; J dt' (fi,Cr) {d i A^ 1 d v C, T ) 
2(3A-1) '* 



(A-l/ 



dt'A- 



2 (a i Ac T )(a l A- 1 a t < T ) + ( + ^ aw ) (^a-^Ct) + (ac t )(sKt) 



A- 1 



tft'A" 



~ {did v h jk ){d% k ) + I (d t >h ij )(Ah i: j) - 3 (didjCrKdeW) 



(B18) 
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We note that the second order gauge transformation is more divergent than the first order one (|B15[) . in the limit 
A — >• 1. 

Finally, we calculate the field £ in the synchronous gauge. Since we adopted a nonperturbative decomposition for 
the spatial metric, it is useful to express this quantity as, 



1 (detg 



Applying the perturbative expansion to the right hand side, we obtain 



12a 2 V" » n cP" y%3 " yv 



Using the transformed metric from (|B12[) . the above expression becomes 



+ C(e 3 ). 



(B19) 



(B20) 



(B21) 



Finally, using the transformations (IB15|) and (|B18p , we obtain Eq. ([57|) , which relates the nonlinear perturbation £ in 
the synchronous gauge to an expansion series of perturbations £t and hij in the transverse gauge. 
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